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ABSTRACT 

We present 10 and 20 /iin images of IRAS 19500-1709 taken with the mid-infrared 
camera, OSCIR, mounted on the Gemini North Telescope. An extended circumsteUar 
envelope is detected, with the N band image indicating an elongation in a NE-SW 
direction. We use a dust radiation transport code to fit the spectral energy distribution 
from UV to sub-mm wavelengths, with a detached dust shell model. A good fit is 
achieved using dust composed of amorphous carbon, silicon carbide and magnesium 
sulphide. We derive estimates for the inner and outer radius, density and mass of the 
dust in the shell. The inner radius is not resolved in our OSCIR imaging, giving an 
upper limit of 0.4 arcsec. With this constraint, we conclude that IRAS 19500-1709 
must be at least 4 kpc away in order to have the minimum luminosity consistent with 
a post-AGB status. 

Key words: radiative transfer - stars: AGB and post-AGB - stars: circumsteUar 
matter - stars: mass loss - stars: individual: IRAS 19500-1709 - stars: individual: HD 
187885. 



1 INTRODUCTION 

Post-AGB stars are evolved stars with low to intermediate 
mass (0.8-8.0 M©) progenitors, that have left the Asymp- 
totic Giant Branch (AGB) but not yet become the central 
stars of Planetary Nebulae (PN) (see review by van Winckel 
2003). The heavy mass- loss that occurred on the AGB has 
ended and the dust envelope that surrounds the star, as a 
result of this mass-loss, expands away from the star, decreas- 
ing in temperature while the temperature of the central star 
increases. During this relatively short (10^-10* years) phase, 
these cool dust rich envelopes radiate most of their energy 
at mid- and far-infrared wavelengths. Observations in the 
mid-infrared probe the inner, warmest part of the envelope, 
the area of most recent mass-loss which occurred near the 
end of the AGB. Studies of individual objects in this transi- 
tional post-AGB phase can be used to investigate this final 
period of AGB mass-loss and to determine the physical and 
chemical properties of the dust. 

IRAS 19500-1709 is associated with the high galac- 
tic latitude F2-3I (Parthasarathy, Pottasch & Wamsteker 
1988) post-AGB star HD 187885 (SAO 163075) and has 
the double-peaked spectral energy distribution (SED) typi- 
cal of post-AGB stars with detached circumsteUar dust en- 
velopes (Hrivnak, Kwok & Volk 1989). The radiation from 
the star, reddened by the dust shell, peaks at around 1 ^m 
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whereas the thermal emission, from the shell itself, peaks 
at 25-30 /im. The post-3rd-dredge-up status of this object 
is confirmed by its high abundance of s-process elements 
(van Winckel & Reyniers 2000). The expansion velocity of 
the envelope, based on the CO-line emission, is 11 km s~^ 
with wings up to 30 km s~^ (Likkel et al. 1987). The enve- 
lope is carbon rich, showing no OH or H2O maser emission 
(Likkel 1989). It has a weak 21 /im feature (Justtanont et al. 
1996) and a broad feature around 30 pm (Hony, Waters & 
Tielens 2002). The origin of these two features is unknown, 
von Helden et al. (2000) identified titanium carbide (TiC) 
as the carrier of the 21/im feature but this has recently been 
disputed by Li (2003). Speck & Hofmeister (2004) have sug- 
gested that SiC may carry the 21 (im. feature. Goebel & 
Moseley (1985) proposed that the 30 /xm feature may be due 
to magnesium sulphide (MgS). The central star has a tem- 
perature of 8000 K (van Winckel & Reyniers 2000). Imaging 
polarimetry has shown that IRAS 19500-1709 has a bipolar 
structure in scattered light (Gledhill et al. 2001). 

In section 2 of this paper we present 10 and 20 /xm im- 
ages of IRAS 19500-1709 taken with the mid-infrared cam- 
era, OSCIR, on Gemini North^. In section 3 we discuss the 



^ Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc; under a cooperative agreement with the NSF 
on behalf of the Gemini partnership:the National Science Foun- 
dation (United States), the Particle Physics and Astronomy Re- 
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results of our modelling of these images and of the SED 
using a dust radiation transport (RT) code. 



2 OBSERVATIONS 

We present 10 and 20 (N and Q3 band) images of IRAS 
19500-1709 taken on 2001 July 12 with OSCIR, mounted on 
the 8.1-m Gemini North Telescope (Fig. 1). OSCIR is a mid- 
infrared camera and spectrometer system built by the Uni- 
versity of Florida Infrared Astrophysics Group. The array 
field of view is 11 x 11 arcsec and the pixel scale 0.089 arcsec. 
We used the standard nodding/chopping mid-infrared ob- 
servational technique to correct for the thermal background 
with a chop frequency of 3 Hz and a chop throw of 15 arcsec. 
For each filter, frame times were 10 ms and on-sourcc expo- 
sure times were 4.3 min. For fiux calibration we used Vega, 
which we observed before and after observing IRAS 19500- 
1709. Atmospheric transmission, as determined by the Vega 
measurements, varied throughout the night and average cal- 
ibration factors have been calculated. The resultant N and 
Q3 band fluxes, obtained by summing over the whole object, 
are given in Table 1. We have also used the Vega observa- 
tions as an indication of the point spread function (PSF) 
and the full width at half-maximum (FWHM) is shown in 
Table 1 along with other observing details. 

Flat-field images were obtained for each filter by expos- 
ing on sky and on a polystyrene flat-field source without 
chopping and nodding and the flat-field constructed using 
the OFLAT task, which is part of the GEMINI package within 
IRAF^. The observations were reduced using the OREDUCE 
task in the GEMINI IRAF package which derives the aver- 
age of the chop and nod differences and produces a single 
fiat-ficld-corrected image for each filter. A few lines of bad 
pixels on these images, caused by an OSCIR channel prob- 
lem which was present during the eaxly Gemini observing 
runs, were fixed using an IDL procedure provided by Scott 
Fisher (OSCIR support team). 

The images of IRAS 19500-1709 show that it is extended 
relative to the PSF, with a FWHM of 1.2 arcsec in both fil- 
ters and diameter of 4 arcsec. The N band image appears 
elongated in the outer contours with position angle (PA) 
of 20 ± 10 dcg. We note that the telescope PSF, as deter- 
mined by observations of Vega, is also elongated in a similar 
direction (shown contoured in Fig. 1). However, the outer 
contours of the Q3 band image appear round, although the 
PSF is also elongated, suggesting that the elongation seen 
in the outer regions of the N band image may be real. This 
interpretation is strengthened by previous lower spatial res- 
olution imaging which also tentatively suggests an elliptical 
structure with PA 30 ± 10 deg. (Meixner et al. 1997). The 
brightness peak in the N band image appears approximately 
0.2 arcsec off centre to the east, relative to the outer con- 
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^ Image Reduction and Analysis Facility is distributed by the Na- 
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Table 1. Measurements for IRAS 19500-1709 for each filter show- 
ing the central wavelength (Ac), bandpass (AA), flux and the 
FWHM of the PSF, estimated from observations of Vega. Due to 
variations in sky transparency, the errors on the fluxes may be as 
large as ±50 per cent. 



Filter 


Ac 


AA 


Flux 


FWHM (Vega) 




(/im) 


(/im) 


(Jy) 


(arcsec) 


N wide 


10.75 


5.23 


33 


0.51 


Q3 


20.80 


1.65 


134 


0.64 



tours. The inner boundary, of any detached shell, is unre- 
solved. 



3 MODELLING AND DISCUSSION 

We use the dust RT code originally developed by Efstathiou 
& Rowan- Robinson (1990), which solves the radiation trans- 
fer equation for an axisymmetric distribution of dust around 
a central source. The code allows up to seven different grain 
materials and has been modified to include a power-law 
size distribution of dust grains (Gledhill & Yates 2003). Al- 
though the code is fully 2-D, we use a spherically symmetric 
dust shell to model 19500-1709 since there is only tenta- 
tive evidence for axisymmetry in the OSCIR observations 
and this would not provide sufficient constraints to justify 
the inclusion of the additional parameters required by a 2-D 
model. 

The input parameters to a spherical shell model include 
the ratio of the stellar radius to the inner shell radius (r*/ri), 
the ratio of the inner shell radius to the outer shell radius 
{ri/r2), and the stellar temperature (Toff). The dust density 
is specified by the extinction (^v) through the shell with the 
dust density assumed to decline with radius as r~^, corre- 
sponding to a constant mass-loss rate. The absorption and 
scattering efficiencies arc calculated for each grain type us- 
ing Mie theory with the dust properties determined by the 
optical constants and grain sizes. We assume spherical dust 
grains with radii between Omin and ttmax with a size distri- 
bution of n{a) DC a~''. 

The resultant model spectrum is plotted against the 
observational data for IRAS 19500-1709, which have been 
taken from the literature. These include optical and near- 
infrared photometry (Hrivnak et al. 1989), a 7.6-23.6 fj,m 
UKIRT spectrum (CGS3) provided by K. Justtanont, IRAS 
fiuxes, 8.5, 10 and 12.2 ^m fiux (Meixner et al. 1997), lUE 
data^ ISO SWS"* and LWS'' spectra and JCMT fiux (van 
der Veen et al. 1994). 

3.1 Fitting the SED 

Around 90 models were run and the best-fitting model, as 
assessed by eye, is shown in Fig. 2. The final model includes 

dust composed of amorphous carbon (amC), silicon carbide 
(SiC) and magnesium sulphide (MgS) grains. The amC dust 
fits the continuum of the SED quite well. 

^ Based on INES data from the lUE satellite 

^ ISO SWS06 observing program rszczerba-PPN30 

^ 750 LWSOl observing program mbarlow-dust 3 
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Figure 1. Images of IRAS 19500-1709 from the OSCIR camera on Gemini Nortii using N and Q3 filters. Tiie PSF is shown contoured 
in the lower right corners. All innermost contours are 95% of the peaJf values and spaced at 0.5 magnitude increments. 



The broad emission feature, extending from 10-13 /im 
and peaking around 12 /xm was observed by Justtanont et 
al. (1996) but has not been conclusively identified. Although 
SiC is thought to be a common dust component in carbon 
stars (Speck et al. 1997), it results in a feature which peaks 
at 11.3 ^m and is too narrow to fit the broad 12 /im 'plateau' 
feature. However, in the absence of an alternative identi- 
fication, we have used a-SiC optical constants (Pegourie 
1988) to produce a feature in this wavelength range. The 8- 
10 lira region is also problematic, in that it exhibits uniden- 
tified features. Kwok, Volk & Bernath (2001) suggest that 
these 8 and 12 /im plateau features may be due to aromatic 
and aliphatic hydrocarbons. Justtanont et al. (1996) pro- 
pose that a poly cyclic aromatic hydrocarbon (PAH), such 
Eis chrysene (C18H12), provides the best match to UIR fea- 
tures in the 8-13 /im region. We have briefiy explored PAH 
materials using the cross-sections of Li & Draine (2001) but 
find that the features produced are too narrow to match the 
plateau-nature of the observed emission and additional fea- 
tures in the 16-25 /iim region of the spectrum are produced, 
which are not seen in the data. 

We include MgS as being the possible carrier of the 
30 pm band. It has been found to provide a reasonable fit 
to this feature (Hony et al. 2002). We used optical constants 
for Mga;Fe(i_a;)S from Begemann et al (1994) with a;=0.9, 
as an approximation to MgS (also see Szczerba et al. 1997, 
Hony et al. 2003). These data do not extend below 10 ^m 
however, so it was necessary to approximate the optical and 
UV extinction cross-sections to ensure that the MgS grains 
absorb sufficient radiation at short wavelengths to radiate in 
the mid-IR. This was done by adopting the optical constants 
of parallel graphite (Draine & Lee 1984) below 4 ^m and 
interpolating between 4 and 10 ^irn to produce a smooth 
extinction curve. The model shows a feature in the UV part 
of the spectrum, between 0.2-0.25 /um. Both parallel graphite 
and amC show structure in their extinction cross sections 
in this region. The feature around 30 /im is broader and 
peaks at a longer wavelength than that produced by our 
model. This may be because we have used spherical grains. 
Variations in grain shape can influence the emission proflle. 



with ellipsoidal grains producing a broader peaJi (Hony et 
al. 2002). Szczerba et al. (1999) suggest that the 30 /xm 
feature may be related to the 21 ^va feature, as stars with 
the 21 /nm feature seem to have the 30 /Ltm feature as well. 
The reverse is not the case though; the 21 //m feature is very 
weak in 19500-1709, which van Winckel & Reyniers (2000) 
suggest may be due to the temperature of this object as the 
feature is observed to be stronger in the cooler 21 pm post- 
AGB stars in their sample. It certainly seems possible that 
the 30 //m feature may be due to an unknown carbonaceous 
component, rather than MgS, as it is a feature which is only 
seen in carbon rich objects. 

In our model the central star is assumed to radiate 

as a black body with Tcff = 8000 K. This temperature is 
higher than has been previously assumed in some of the 
literature (e.g. Hrivnak et al. (1989), Giirtler, Kompe & 
Henning (1996), Meixner et al. (1997) use temperatures 
of 7000-7500 K which have been estimated from the spec- 
tral type). However detailed analysis of spectral line data 
suggests a higher temperature of 8000 K (van Winckel & 
Reyniers (2000)) so we adopt this as the most reliable es- 
timate. Paxthasarathy et al. (1988) suggested a tempera^ 
ture of 8500 K might be required to fit the UV spectrum of 
IRAS 19500-1709 which they found to be peculiar because 
of a broad absorption feature around 0.17 /im. Increasing 
the temperature of the central star improves the fit around 
the UV part of the SED but also increases the prominence 
of the SiC emission feature. So an increase in temperature 
needs to be balanced with a decrease in the abundance of 
SiC. However if the abundance of amC is too high the model 
does not fit the data around 15-25 /um. A dust mixture com- 
prising 40 per cent amC, 30 per cent SiC and 30 per cent 
MgS provided the best fit to the data. 

A value of Av=1.2 provided the best fit around the UV 
and optical part of the spectrum. Increasing A^, results in 
too much extinction in this region and also too much mid- 
IR flux. Larger dust grains produce too much flux in the 
UV and optical region and too strong an emission feature 
around 11 /xm but insufficient flux in the mid-IR peak. A 
value of Omin=0.01 /xm yielded the best results along with 
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Figure 2. SED and model fit (solid black line). 

a grain size power-law index of g=5. Reducing the value of 
q includes too many larger grains and so also results in a 
mid-IR peak which is too weak. 

The fit to the 10-25 /im region is strongly affected by 
r«/ri which determines the temperature of the warmest 
dust. A value of 8 x 10~^ provided the best fit. Increas- 
ing this value results in dust grains which are too hot and 
a poor fit around this region. Reducing the value of rt/ri 
results in dust which is too cold and so not enough flux at 
mid-IR wavelengths and too much in the 40-160 fjia region. 
The ratio of the inner shell radius to the outer shell radius 
{ri/ri) is determined to be 0.1 as this produces sufficient 
far-IR emission to fit the data longward of 100 iim.. Extend- 
ing the boundary any further than this produces too much 
far-IR emission. 

After fitting the model to the data, the model curve was 
scaled to fit the observed fiux. Table 2 shows the parameters 
for the model fit along with the derived parameters. 

3.2 Derived parameters and the distance to 
19500-1709 

The profiles of the raw model images are double-peaJsed 
showing a resolved inner boundary. The inner boundary 
of the OSCIR images, however, is unresolved. In order to 
match the model to the data and so determine the scale 
of the model (i.e. the size of a model pixel in arcsec), the 
model images were smoothed using a Gaussian filter until 
they resembled the OSCIR images. Using smoothing steps 
of 1 model pixel, the Gaussian smoothing PSF was deter- 
mined to have FWHM > 11 pixels. A value less than this 
resulted in a dip in the profile so that FWHM=11 is the low- 
est value to show a single peak in the profile. Equating this 



degree of smoothing to the observed PSF (see Table 1) gives 
a model pixel size of 0.046 arcsec. This then leads to values 
of ri < 0.4 arcsec and r2 < 4 arcsec (for ri/r2 = 0.1). Since 
r,/ri = 8 X 10"^ then r, < 6.87 R© and L < 173 Lq. These 
and other resulting upper limits on the derived parameters 
are listed in Table 2 for a distance of 1 kpc. 

Schonberner (1983) determined lower limits for the 
mass and luminosity of central stars of PN of approximately 
0.55 Mq and 2500 L©. If IRAS 19500-1709 is to be this lumi- 
nous then it must be at a distance of at least 4 kpc. Pottasch 
& Parthasarathy (1988) conluded that D = 5.5 kpc based 
on the optical-IR energy balance. The value of r,/ri is well 
constrained by the model fit: increasing or decreasing this 
value results in dust that is either too hot or too cool, re- 
spectively. We explored the possiblility that axisymmetric 
models might allow us to vary r*/ri, by redistributing the 
dust into a smaller volume of the shell, however the effect 
was minimal. Meixner et al. (1997) have used a larger value 
of r,/ri = 2.7 X 10"" but it is evident from their SED fit 
that the dust is too hot. With the model constraining r*/ri 
in this way we conclude that IRAS 19500-1709 must be at 
least 4 kpc away for it to have the minimum luminosity con- 
sistent with a post-AGB status and this distance is adopted 
in the rest of the paper. 

Dividing ri by the expansion velocity of the envelope 
(11 km s"\ Likkel et al. 1987), shows that, for D = 4 kpc, 
the shell became detached (i.e. the high mass loss phase 
ended) at least 690 years ago. This is comparable to the 
result of 740 years obtained by Meixner et al. (1997). The 
age of the shell, calculated by dividing r2 by the expansion 
velocity, is at least 6900 years. 

The mass of dust in the shell and the grain number den- 
sity at ri were calculated using an average of the densities 
of the 3 grain types used in the model and an average value 
for the extinction cross sections. Using the dust shell dimen- 
sions derived from the model fit, then to produce an extinc- 
tion through the dust shell of = 1.2, which is needed 
to fit the SED, requires a dust number density at the inner 
shell radius of ni = 5.6 x 10~^ cm~^, cissuming f = 4 kpc. 
Using this number density, assuming a dust density dis- 
tribution and an average grain bulk density of 2.56 g cm~^ 
results in a total mass of dust of Md = 1.89 X 10"^ Mq. 
Gledhill, Bains & Yates (2002) used their 850 ^m SCUBA 
photometry to obtain a lower limit of Md > 1.6 x 10~^ M© 
at 4 kpc. 

Using a dust-to-gas mass ratio of 4.5 x 10"'^, an approx- 
imate value used for carbon-rich AGB stars (Jura 1986), the 
total envelope mass at D = 4 kpc is 0.42 M© . According to 
our calculations the mass loss episode, responsible for the 
shell, lasted 6200 years so this results in an average mass 
loss rate of 6.8 x 10"'' M© yr"^ at this distance. This is 
comparable to the mass loss rate (Al) derived by Giirtler et 
al. (1996) of 6 X 10"^ M© yr"^ using a RT code and amC 
dust grains. It lies between the values derived by Meixner 
et al. (1997, M = 1.4 x 10"" M© yr"^) and Hrivnak et al. 
(1989, M = 1.5 X 10"'^ M© yr'^). Omont et al. (1993) cal- 
culate a value of M = 1.2 - 1.5 x 10"^ M© yr"^ from CO 
line emissions. All of these rates are for D = 4 kpc. 
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Table 2. Model and derived parameters resulting from the final 
fit to the SED. 



Parameter 



Value 



r*/ri 

ri/r2 
Av (mag) 

q 

Imin (^in) 

(^tm) 
Teff (K) 
Grain type 
amC (%) 
SiC (%) 
MgS (%) 
Grain temp, (min) 
amC (K) 
SiC (K) 
MgS (K) 
Grain temp, (max) 
amC (K) 
SiC (K) 
MgS (K) 
Derived parameters 
ri (cm/kpc) 
r2 (cm/kpc) 
r, (Ro/kpc) 
L (Lc5/kpc2) 
ni (cm~'^/kpc~^) 
Md (M0/kpc2) 



8x10"° 
0.1 
1.2 

5 
0.01 
2 

8000 

40 
30 
30 

53 
48 
42 

170 
155 
125 

< 5.98x1016 

< 5.98x101*^ 

< 6.87 

< 173 

> 2.2x10""' 

< 1.16x10-4 



3.3 Evidence for axisymmetry 

Ncar-IR imaging polarimctry observations by Glodhill ct al. 
(2001) show evidence for two scattering peaks on either side 
of the star, which suggest scattering from the inner edge 
of an axisymmetric envelope. The hne joining the peaks is 
at PA~ 100°, approximately orthogonal to the elongation 
in the N band image of Fig. 1. The morphology appears 
remarkably similar to that described in mid-IR and near- 
IR imaging of 17436+5003 (Gledhill & Yates 2003; Gledhill 
et al. 2001). Although there is insufficient information in 
our OSCIR images to justify axisymmetric RT models, we 
suggest that 19500-1709 has an optically thin detached ax- 
isymmetric dust shell, resulting from enhanced equatorial 
mass-loss. 



4 CONCLUSIONS 

We present mid-infrared images and a model of the SED of 
IRAS 19500-1709 using amC, SiC and MgS grains. We find 
that small dust grains (amin = 0.01 jim) and a steep power 
law size distribution (^=5) provide the best fit to the data. 
The observations constrain the inner radius of any detached 
dust shell to be a maximum size of 0.4 arcsec. Our results 
suggest that this object must be at a distance of at least 4 
kpc to be sufficiently luminous to be a post-AGB star. We 
show that, at this distance, IRAS 19500-1709 would have an 
envelope mass of 0.42 M©, that during the high mass-loss 
phase the mass-loss rate was 6.8 x 10~^ Mq yr"! and that 
this phase ended ~ 700 yrs ago. Although the inner enve- 
lope structure has not been resolved in these OSCIR ob- 



servations, there is evidence for axisymmetry along an axis 
perpendicular to that seen in near-IR imaging observations. 
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